# -*- coding: utf-8 -*-
"""
Spyder Editor

This is a temporary script file.
"""
from NewmarkCD import newmarkCD 
import numpy as np

def spectrumD(file,dt,Npt,T):
    
    fa = np.loadtxt(file)*9805
    
    Sd = np.zeros(len(T))
    
    for i in range(1,len(T)):
        Sd[i] = newmarkCD(T[i],0.05,fa,dt,Npt)/1000
        
    return Sd
    




